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ABSTRACT 


Describes  results  of  testing  a  generalized  equation  for  predicting 
bark  thickness  using  data  from  western  Montana  species.     The  most 
satisfactory  least  squares  solution  for  use  with  a  combination  of 
species  is  given. 


Upper-stem  diameters  measured  using  an  optical  dendrometer  include  the  bark. 
Therefore,  to  estimate  volume  of  wood  alone,  the  measurement  must  be  adjusted  to  diameter 
inside  bark.    However,  average  bark  thickness  and  the  pattern  of  its  distribution  on  the  stem 
may  vary  among  species.    It  might  be  possible  to  develop  a  satisfactory  equation  for  predicting 
bark  thickness  for  each  species  in  an  area  if  a  large  enough  sample  of  bark  thickness  measure- 
ments were  available  for  each  species.   When  such  information  is  not  available,  it  would 
seem  reasonable  to  use  a  generalized  model  for  predicting  bark  thickness  in  which  one  of  the 
variables  is  a  measurement  of  bark  thickness  at  breast  height  on  the  standing  tree.  This 
variable  should  explain  much  of  the  difference  in  upper-stem  bark  thickness  between  species 
and  between  trees,  because  the  factors  that  affect  bark  thickness  (such  as  age  and  species)  up 
the  stem  also  affect  thickness  at  breast  height. 

Grosenbaugh  (1967)  incorporated  three  options  for  estimating  bark  thickness  in  his  STX 
program  for  converting  dendrometer  measurements  on  standing  trees  to  volume,  surface  area, 
and  length.    One  of  these  was  the  equation 
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bark  thickness  at  the  /'   level  on  the  tree  stem 
diameter  outside  bark  at  the  i*   level  on  the  tree  stem 
bark  thickness  at  breast  height 
diameter  at  breast  height  outside  bark 


This  equation  can  also  be  used  to  estimate  bark  thickness  in  the  NETVSL  computer  program, 
which  was  designed  to  calculate  volumes  of  felled  or  standing  trees  (Stage,  Dodge,  and  Brickell 
1968).   Whether  single  or  double  bark  thickness  is  used  is  immaterial.    If  B4,5is  a  single  bark 
thickness  measurement,  the  predictions,  Bi5  will  be  for  single  bark  thickness.    If  B4-6is  a 
double  bark  thickness  measurement,  then  predictions  will  be  of  double  bark  thickness.  All 
fractions  in  the  equation  are  equal  to  unity  and  predicted  bark  thickness  is  equal  to  actual  bark 
thickness  when  Dj  is  D4  B. 

Mesavage  (1969)  found  that  equation  (1)  gave  reasonably  accurate  bark  thickness  predic- 
tions for  four  southern  pine  species  of  sawtimber  age.    For  other  species  and  ages,  where 
predictions  from  equation  (1)  may  not  be  as  accurate,  we  can  use  the  equation 


where  Q  is  a  coefficient  to  be  estimated  from  a  sample  of  bark  thickness  measurements  taken 
from  the  population  for  which  the  equation  will  be  used.    The  coefficient  Q  can  be  estimated  by 
least  squares  regression. 

Even  more  flexibility — perhaps  too  much — will  result  from  including  more  coefficients 
so  that 


where  the  b; 's  are  coefficients  to  be  estimated  and  bg  replaces  Q  in  equation  (2).    Equation  (3) 
retains  the  property  of  predicting  a  bark  thickness  equal  to  actual  bark  thickness  when  D;  is 
D4.6-    Equation  (2)  is  equivalent  to  equation  (3)  when  coefficients  b±,  b3,  and  b4  are  restricted 
to  unity.    In  this  paper,  we  will  describe  the  results  of  testing  equations  (2)  and  (3). 


Using  trees  felled  on  the  Kootenai  National  Forest  in  western  Montana,  we  took  443 
measurements  of  double  bark  thickness.    The  following  species  were  represented  in  the  sample 
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DATA  USED  FOR  TESTING  THE  MODEL 


Western  white  pine 
Ponderosa  pine 
Lodgepole  pine 
Grand  fir 
Subalpine  fir 
Inland  Douglas-fir 


(Pinus  monticola  Dougl. ) 

(Pinus  ponderosa  Laws. ) 

(Pinus  contorta  var.  latifolia  S.  Wats.) 

(Abies  grandis  Lindl. ) 

(Abies  lasiocarpa  Nutt. ) 

(Pseudotsuga  menziesii  var.  glauca 


Engolmann  spruce 
Western  hemlock 
Western  redcedar 


Western  larch 


(Mirb. )  Franco) 
(Larix  occidentalis  Nutt. ) 
( Picea  cngclmannii  Parry) 
(Tsuga  heterophylla  (Raf. )  Sarg. ) 
(Thuja  plicata  D.  Don) 
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Measurements  were  taken  on  sawn  cross-sections  of  the  trunk  at  points  8  feet  apart,  beginning 
4.6  feet  above  breast  height  and  continuing  up  the  stem  until  a  point  was  reached  where  diameter 
outside  bark  was  less  than  4  inches.    All  bark  measurements  used  were  taken  from  points 
where  diameter  exceeded  4  inches.    The  maximum  and  minimum  axes  of  the  cross-section  were 
measured  outside  bark  and  called  DOB1  and  DOB2,  respectively.    Inside  bark  measurements  on 
the  same  axes  were  called  DII^  and  DIB2.   Then  double  bark  thickness  was  calculated  as 

DBT  =  VDOB^OBg   -  V  DIB1DIB2  . 

All  measurements  were  taken  to  the  nearest  1/100  foot  (0.12  inch),  corresponding  to  single 
bark  thickness  measurement  to  the  nearest  0.06  inch.    The  average  double  bark  thickness  over 
the  entire  sample  was  0.  96  inch. 

USE  OF  COMPUTER  PROGRAM  "NLLS"  TO  FIT  NONLINEAR  REGRESSIONS 

Equations  (2)  and  (3)  cannot  be  made  linear  with  respect  to  their  parameters  by  transfor- 
mation.   Therefore,  they  cannot  be  fitted  to  data  using  linear  regression  methods.  Instead, 
the  unknown  coefficients  were  estimated  using  a  computer  program  for  nonlinear  regression 
(Marquardt  1963,  1966;  Draper  and  Smith  1966)  that  minimizes  the  sum  of  squared  residuals. 
To  use  the  NLLS  (nonlinear  least  squares)  program,  the  desired  function  is  written  into  a 
FORTRAN  subroutine.    Initial  estimates  must  be  provided  for  each  coefficient  to  start  the 
iterative  fitting  process. 

FITTING  THE  BARK  THICKNESS  MODEL 

Equation  (2)  was  fitted  first.    This  is  a  one-parameter  model  and  the  initial  estimate  of  the 
coefficient  was  2.0.    This  made  the  equation  equivalent  to  equation  (1)  and  made  it  possible  to 
see  what  improvement  in  fit,  if  any,  could  be  obtained  by  a  different  estimate  of  the  coefficient. 
The  result  was  as  follows: 

Starting  Final  solution 

2.0  1.67342 
86.8995  80.2518 
.4434  .4261 

The  reduction  in  sum  of  squares,  while  statistically  significant  according  to  the  test  described 
by  Williams  (1959,  p.  81),  is  small.    The  standard  error  of  residuals  was  reduced  by  less  than 
0.  02  inch. 

Because  it  is  more  flexible,  equation  (3)  could  be  expected  to  account  for  more  variation 
than  was  possible  with  equation  (2).    The  results  of  fitting  equation  (3)  were: 


Starting  Final  solution 


Coefficient  estimate  bx 

ho 

0.548235 

Coefficient  estimate  b2 

2.0 

5. 80698 

Coefficient  estimate  b3 

1.0 

4. 84516 

Coefficient  estimate  b4 

1.0 

4.40345 

Sum  of  squared  residuals 

86.  8995 

68.  7208 

Standard  error  of  residuals 

.4449 

.  3957 

Coefficient  estimate 
Sum  of  squared  residuals 
Standard  error  of  residuals 
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These  results  show  a  more  substantial  reduction  in  residual  sum  of  squares  than  was  obtained 
with  equation  (2).    However,  the  high  correlation  between  the  estimates  of  the  parameters  indi- 
cated that  a  two-parameter  model  might  provide  an  equally  close  fit.    Therefore,  coefficients 
b3  and  b4  were  deleted  from  equation  (3);  that  is,  they  were  held  to  a  constant  value  of  unity, 


giving 


Fitting  this  equation  gave: 


B;    =  B< 


(d4.b) 


bP  -  1 


bP~ 


(  D4.5) 
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Final  solution 


Coefficient  estimate  bx  0.0352 
Coefficient  estimate  bz  1.  25638 

Sum  of  squared  residuals  69.4400 
Standard  error  of  residuals  .  396813 


Using  this  model,  the  coefficient  of  correlation  between  estimates  of  the  coefficients  was 
0.  8958. 


Variance  around  the  surface  described  by  equations  (3)  and  (4)  was  very  nearly  propor- 
tional to  observed  bark  thickness.    Equation  (4)  was  fitted  a  second  time  assigning  each  obser- 
vation a  weight  inversely  proportional  to  bark  thickness  so  as  to  attain  the  condition  of  homo- 
gencouH  variunco  about  regression.   The  resulting  estimates  were: 


Final  solution 


-0. 009363 

1. 22228 
70. 9300 
.401047 

Two  observations  can  be  made  from  these  figures.    First,  weighting  of  the  observations 
did  not  result  in  much  change,  either  in  the  estimated  coefficients  or  in  the  residual  sum  of 
squares.    In  other  words,  weighting  did  not  change  the  position  or  shape  of" the  regression  sur- 
face appreciably.   This  may  indicate  that  the  equation  being  fitted  is  quite  appropriate  to  the 
relationship  between  variables  represented  by  these  data.    Second,  it  is  noteworthy  that  the 
coefficient  bx  is  nearly  zero.    This  indicates  the  value  for  b^  that  reduces  the  variation  in  the 
quantity 


Coefficient  estimate  bx 
Coefficient  estimate  bP 

2/ 

Sum  of  squared  residuals-' 
Standard  error  of  residuals 


to  a  minimum  best  satisfies  the  requirements  of  least  squares.    The  above  expression  would 
always  be  equal  to  one,  if  bx  were  exactly  zero,  no  matter  what  the  values  of  Dj  and  D4  5  .. 
Variation  in  the  ratio  of  diameter  outside  bark  to  diameter  at  breast  height  (as  the  ratio  appears 
with  an  exponent  in  equation  (4) ),  apparently  does  not  explain  any  variation  in  observed  bark 
thickness.    Rather,  variation  in  the  D;  /D4.6  ratio  seems  to  reduce  precision  of  the  predicting 

g  Unweighted  sum  of  squares  and  standard  error  of  deviations  from  the  regression  surface 
fitted  to  weighted  observations.    The  figures  are  presented  in  this  way  so  as  to  be  comparable  to 
those  from  unweighted  regressions. 
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equation  except  for  its  use  in  the  denominator  of  the  fraction.  If  bx  is  assumed  equal  to  zero, 
then  the  ratio  raised  to  that  power  will  always  be  equal  to  one,  and  equation  (4)  becomes 


Bj    -  B4J 


b2  -  1 


D; 


Equation  (5)  is  still  a  special  case  of  equation  (3). 
inversely  proportional  to  bark  thickness  it  gave: 


(5) 


Fitted  with  observation  weights 


Final  solution 

Coefficient  estimate  b2  1.  23177 

Sum  of  squared  residuals  70. 1327 


Predictions  of  double  bark  thickness  given  by  this  equation  are  shown  in  figure  1.    Double  bark 
thickness  would  logically  be  expected  to  approach  zero  as  the  ratio  of  upper-stem  diameter  to 
diameter  at  breast  height  approaches  zero.    Equation  (5)  does  not  meet  this  requirement.  None 
of  the  bark  thickness  measurements  used  to  derive  the  equation  were  taken  from  points  on  the 
tree  stem  where  diameter  outside  bark  was  less  than  4  inches.    If  bark  measurements  associ- 
ated with  smaller  diameters  had  been  used  it  seems  quite  possible  that  a  properly  fitting 
equation  would  have  been  required  to  follow  the  marked  decrease  in  bark  thickness  that  usually 
occurs  in  the  top  of  the  tree.    For  practical  use  this  seems  to  be  rather  unimportant,  because 
material  In  treelopB  wmallcr  than  4  inches  in  diameter  is  usually  unmerchantable. 


DIAMETER  INSIDE  BARK/DIAMETER  OUTSIDE  BARK  RATIO 


Mesavage  (1969)  cited  other  authors'  findings:  that  this  ratio  was  nearly  constant  at  all 
levels  on  the  stem  for  some  species;  that  it  increased  with  increasing  height  for  other  species; 
and  that  it  decreased  with  increasing  height  for  still  other  species.    MacDonald  (1933)  examined 
the  trend  of  the  DIB/DOB  ratio  up  the  stem  for  several  coniferous  species  in  Great  Britain.  He 
found  that,  in  general,  the  ratio  decreases  with  increasing  distance  above  breast  height  through 
the  lower  third  of  the  tree  bole,  remains  nearly  constant,  then  increases  in  the  upper  half  of 
the  stem.    According  to  the  bark  distribution  pattern  described  by  equation  (5),  the  DIB/DOB 
ratio  increases,  remains  nearly  constant,  then  decreases  as  measurement  progresses  to 
points  further  up  the  tree  stem.     The  DIB/DOB  ratios  at  various  points  on  the  stem  of  a  tree 
having  a  diameter  of  10  inches  and  a  double  bark  thickness  of  1  inch  at  breast  height  follow: 


Diameter 

Double  bark  thickness 

DIB/DOB 

103 

31.00 

0.  900 

9 

.70 

.  922 

8 

•  .54 

.  933 

7 

.44 

.938 

6 

.37 

.  939 

5 

..32 

.  937 

4 

.28 

.  930 

3 

.25 

.  917 

2  ■ 

.22 

.  888 

1 

.20 

.  795 

3  At  breast  height. 
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RATIO  OF  DIAMETER  OUTSIDE  BARK  TO  DIAMETER  AT  BREAST  HEIGHT 

(D  (I)  /D  (4.5) ) 


Figure  1. — Double  bark  thickness  predicted  from  given  measurements  of  double  bark  thick- 
ness at  breast  height  and  ratios  of  diameter  outside  bark  to  diameter  at  breast 
height. 


However,  if  the  bark  thickness  model  given  by  equations  (2)  and  (3)  is  appropriate  and 
will  accommodate  either  the  increasing,  decreasing,  or  constant  trend  in  DIB/DOB  ratio, 
then  fitting  either  equation  to  bark  thickness  measurements  by  least  squares  will  give  estimates 
of  the  coefficients.    A  great  deal  of  concern  over  the  DIB/DOB  ratio  itself  does  not  seem  . 
justified  if  an  equation  is  available  that  provides  satisfactory  predictions  of  bark  thickness  for 
individual  species  or  groups  of  species.  It  is  bark  thickness,  not  the  ratio,  that  is  the  item  of 
ultimate  interest. 

CONCLUSIONS 

Adequate  information  might  allow  derivation  of  an  equation  to  predict  bark  thickness  for 
each  species.    This  could  lead  to  greater  accuracy  and  precision  in  prediction  if  the  distribu- 
tion of  bark  on  the  tree  stem  varies  from  one  species  to  another.    Possibly,  the  ratio  of  diam- 
eters left  out  of  equation  (5)  should  be  included  for  some  species  as  it  was  in  equations  (2)  and 
(4).   For  some  species,  the  coefficient  values  given  in  equation  (1)  might  give  the  best  fit  by  the 
least  squares  criterion.    Even  if  equation  (1)  were  most  appropriate  for  some  particular  species, 
it  appears  that  equation  (5)  would  be  superior  for  general  application  to  all  species  in  the 
vicinity  of  the  Kootenai  National  Forest.    One  would  hope  that  in  a  population  to  which  the 
equation  might  be  applied,  each  species  would  be  represented  in  the  same  proportion  as  in  the 
sample  on  which  the  equation  was  based.    The  argument  that  this  might  not  be  the  case  mitigates 
against  any  other  all-species  equations  fully  as  much  as  against  equation  (5). 
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